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LONG-TERM  GOAL 

To  develop  a  methodology  for  remotely  measuring  the  spatial  and  temporal  variability  of  the 
atmospheric  surface  latent  and  sensible  heat  fluxes  over  the  Labrador  Sea  using  passive  microwave 
imagery  (SSM/I). 

OBJECTIVES 

(1)  To  test  existing  algorithms  for  calculating  integrated  water  vapor  (IWV)  and  surface  wind  speed 
and  compare  these  calculations  to  shipboard  measurements  made  from  the  Knorr  during  the  1997 
experiment.  Calculate  surface  flux  fields  over  the  Labrador  Sea. 

(2)  Develop  algorithms  optimized  for  use  in  the  Labrador  Sea  using  the  1997  data  and  apply  these  to 
the  1998  data  set. 

APPROACH 

It  has  been  well  proven  that  one  can  estimate  atmospheric  integrated  water  vapor  (IWV)  and  surface 
wind  speed  using  SSM/I  data.  Algorithms  have  been  developed  by  using  oceanic  sounding  data  sets 
and  buoys  in  conjunction  with  radiative  transfer  models  and  multiple  regression  techniques  to  relate 
brightness  temperatures  (Tb’s)  from  different  SSM/I  channels  and  IWV  and  wind  speed  at  19.5  m 
height.  The  SSM/I  instrument  on  the  DMSP  series  of  satellites  operates  at  4  frequencies  (19,  22,  37 
and  85  GHz),  and  at  horizontal  (H)  and  vertical  (V)  polarizations  at  19,  37  and  85  GHz  and  at  V 
polarization  at  22  GHz.  Spatial  resolution  is  nominally  50  km  at  19,  22  and  37  GHz  and  25  km  at  85 
GHz. 

Using  over-the-ocean  radiosondes  Liu  (1986)  showed  that  IWV  was  highly  correlated  with  surface 
mixing  ratio  (qa)  on  a  monthly  averaged  basis.  Subsequently,  Hsu  and  Blanchard  (1988)  showed  this 
relationship  held  for  individual  soundings.  Previous  studies  have  shown  that  latent  heat  fluxes  can  be 
estimated  from  quantities  (IWV  and  wind  speed)  that  can  be  estimated  from  SSM/I  data  (Claud  et  al, 
1992;  Clayson  and  Curry,  1996). 

Claud  et  al  (1992)  developed  an  empirical  relationship  between  SSM/I  brightness  temperatures,  an 
independent  estimate  of  Tsfc,  and  (qs-qa),  where  qs  is  the  saturation  mixing  ratio  at  the  ocean  surface 
temperature  Ts^,  qa  is  the  mixing  ratio  of  the  air  near  the  ocean  surface.  Having  (qs-qa)  and  wind  speed 
allows  us  to  calculate  a  value  of  the  latent  heat  flux,  LHF,  using  the  bulk  method: 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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WORK  COMPLETED 

I  followed  the  above  approach  in  developing  a  technique  for  estimating  LHF  over  the  Labrador  Sea 
during  the  1997  field  program.  The  212  radiosondes  launched  from  the  Knorr  during  February  2  - 
March  13,  1997  provided  the  in  situ  data  from  which  I  developed  the  algorithms  for  calculating  IWV 
and  wind  speed  and  ultimately  LHF. 

I  developed  the  relationship  between  IWV  and  T22v  and  T19v  by  regressing  the  values  of  the 
integrated  water  vapor  calculated  from  the  Knorr  soundings  and  the  brightness  temperature  values  of 
T22v  and  T19v  at  the  Knorr  location.  All  the  following  relationships  were  developed  on  a  daily- 
averaged  basis.  The  resulting  IWV  algorithm  is: 

IWV  (g/kg)  =  -16.4415  -I-  .9302*T22v  -  .8563*TI9v. 

I  followed  the  form  of  the  wind  speed  algorithm  developed  by  Goodberlet  et  al  (1989)  where 

u  =  ao  -t  ai*T19v  -i-  a2*T22v  -i-  a3*T37v  -i-  a4*T37h, 

where  the  constants  ao  -  a4  are  obtained  from  a  regression  analysis. 

Using  the  values  of  the  wind  speed  adjusted  to  a  height  of  10  m  measured  on  the  Knorr  and  the  SSM/I 
values  of  T19v,  T22v,  T37v  and  T37h  at  the  location  of  the  Knorr  I  developed  the  following  wind 
speed  algorithm: 

u  (m/s)  =  193.557  -i-  .3777*T19v  -  .0276*T22v  -i-  L7815*T37v  -i-  .8195*T37h 

Following  the  approach  of  Claud  et  al  (1992),  using  the  calculated  value  of  qa(at  10  m)  from  the 
sounding  data  and  qs  calculated  from  the  measured  values  of  Tscf  at  the  Knorr  I  similarly  developed  the 
following  algorithm  for  (qs-qa): 

(qs-qa)  =  3.213  -  .2605*IWV  -t  .33104*T,fe. 

In  summary  the  procedure  to  calculate  LHF  is  as  follows: 

(1)  IWV  =  f(T22v,  T19v) 

(2)  u  =  f(TI9v,  T22v,  T37v,  T37h) 

(3)  (qs-qa)  =  f(IWV,  T.^) 

(4)  LHF  =  Lv*p*CE*u*(qs-qa),  (in  W/m^) 

SSM/I  brightness  temperature  fields  were  obtained  for  the  area  5 1°  N  -  67  °  N  and  67  °  W  -  43  °  W, 
containing  the  Labrador  Sea.  Using  the  above  algorithms  for  IWV,  (qs-qa)  and  u,  I  then  calculated  LHF 
fields  for  this  area. 

An  estimate  of  the  sensible  heat  flux,  SHF,  can  also  be  made  from  SSM/I  data  using  the  fact  that  the 
air  temperature  Ta  is  related  to  the  value  of  the  integrated  water  vapor  (Jourdan  and  Gautier,  1995). 


Using  the  daily  averaged  values  of  Ta  from  the  Knorr  and  IWV  from  the  radiosondes  I  obtained  the 
following  relationship  between  Ta  and  IWV : 


Ta  =  -25.27*exp(-.305*IWV)  +  2.0. 

Using  the  bulk  flux  method  the  SHF  can  then  be  calculated  from: 
SHF  =  p*Cp*CH*u*(T,fe-Ta),  (in  W/m"). 


RESULTS 

A  methodology  has  been  developed  for  estimating  latent  and  sensible  heat  flux  fields  over  the 
Labrador  Sea  using  remote  sensing  passive  microwave  data  from  SSM/I.  This  allows  the  specification 
of  spatial  and  temporal  changes  in  the  surface  heat  budget  over  the  wintertime  Labrador  Sea  at  spatial 
scales  (~  50  km)  much  higher  than  is  possible  with  NCEP  reanalysis  fields.  Having  these  higher 
resolution  flux  fields  is  critical  for  studying  the  scales  of  oceanographic  processes  important  to  deep 
convection  in  the  Labrador  Sea. 

Latent  heat  flux  fields  have  been  calculated  for  the  Labrador  Sea  and  qualitatively  show  spatial 
variations  that  one  would  expect.  This  is  also  true  of  the  sensible  heat  flux  fields.  The  average  Bowen 
ratio  (SHF/LHF)  for  the  measurements  on  the  Knorr  was  1.02  whereas  the  Bowen  ratio  for  the 
predicted  fluxes  was  1.24. 

The  correlation  coefficient  between  latent  heat  fluxes  obtained  from  the  in  situ  measurements  on  the 
Knorr  and  those  obtained  from  the  SSM/I  algorithm  was  r  =  0.907.  The  correlation  coefficient 
between  the  Knorr’s  sensible  heat  flux  and  that  from  the  SSM/I  algorithm  was  r  =  0.88. 

The  latent  and  sensible  heat  flux  fields  have  been  obtained  from  the  NCEP  reanalyses  for  the  same 
time  period  (Peb.2  -  Mar.  13)  and  will  be  compared  to  the  flux  fields  estimated  from  the  SSM/I  data. 
The  fluxes  will  also  be  calculated  for  the  period  and  location  of  the  Knorr  during  the  1998  cruise  and 
will  be  compared  to  the  bulk  estimates  made  from  the  in  situ  ship  measurements.  Spatial  fields  of 
fluxes  will  also  be  calculated  for  1998  and  these  fields  compared  to  the  NCEP  reanalysis  flux  fields  for 
the  same  time  period. 

IMPACT/APPLICATION 

If  this  method  of  remotely  estimating  surface  turbulent  sensible  and  latent  heat  flux  fields  proves 
robust  it  would  provide  a  method  for  supplying  data  for  initializing  oceanographic  numerical  models, 
and  for  the  analysis/interpretation  of  oceanographic  measurements.  Having  these  higher  resolution 
(~50  km)  flux  fields  is  critical  for  studying  the  scales  of  oceanographic  processes  important  for  deep 
convection  in  the  Eabrador  Sea. 

TRANSITIONS 

The  above  technique  can  be  applied  to  other  geographic  areas  for  similar  wintertime  cold  air  outbreak 
conditions. 


RELATED  PROJECTS 


(1)  I  have  been  working  with  Kent  Moore  at  the  University  of  Toronto  in  the  interpretation  and 
analysis  of  the  NCEP  reanalysis  fields.  He  is  also  supplying  me  with  AVHRR  imagery  over  the 
Labrador  Sea  recorded  during  the  1997  field  program. 

(2)  I  will  be  comparing  the  heat  flux  fields  that  I  obtain  with  those  that  Russ  Davis  at  SIO  obtained 
from  his  drifter  measurements. 

(3)  I  will  also  be  comparing  my  fluxes  with  those  Eric  D’Asaro  at  APE,  University  of  Washington 
obtained  from  his  Lagrangian  floats. 
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